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Abstract 

The  so-called  brick  layer  model  is  frequently  used  to  analyze  impedance  spectra  of  polycrystalline  samples  with  highly 
resistive  grain  boundaries.  However,  the  basic  assumptions  of  the  model  (cubic  grains,  laterally  homogeneous  grain 
boundaries,  identical  properties  of  all  grain  boundaries)  are  usually  violated  in  real  ceramics.  To  investigate  the  impact  of 
some  deviations  from  the  brick  layer  model,  the  potential  distributions,  and  thus  the  impedance  of  polycrystals,  have  been 
calculated  by  the  finite  element  method.  The  results  show  that  bulk  properties  can  distinctly  influence  the  size  and  shape  of 
the  so-called  ‘grain  boundary  semicircle’,  particularly  for  laterally  inhomogeneous  grain  boundaries  and  for  properties 
varying  from  boundary  to  boundary.  Depressions  of  the  grain  boundary  semicircle  solely  due  to  a non-brick-layer 
microstructure  are  observed.  The  validity  and  limits  of  the  brick  layer  model  are  discussed.  © 2000  Elsevier  Science  B.V. 
All  rights  reserved. 
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1.  Introduction 

Highly  resistive  grain  boundaries  are  of  great 
technological  importance  with  regard  to  the  overall 
electrical  properties  of  polycrystalline  materials  and 
can  play  a desired  functional  role  in  semiconducting 
electroceramic  devices  such  as  varistors,  thermistors 
and  capacitors  [1,2].  On  the  other  hand,  for  many 
applications  of  solid  electrolytes  (e.g.  fuel  cells, 
oxygen  pumps  or  electrochemical  reactors)  the  over- 
all conductivity  of  the  electroceramic  material  should 
be  as  high  as  possible  and  highly  resistive  grain 
boundaries,  such  as  those  found  in  oxygen  conduct- 
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ing  membranes  of  zirconia  [3-7],  are  distinctly 
unwelcome. 

Impedance  spectroscopy  is  an  important  tool  to 
investigate  the  electrical  properties  of  highly  resistive 
grain  boundaries  since  the  grain  boundaries  frequent- 
ly cause  an  additional  semicircle  in  the  complex 
impedance  plane  [3-7],  In  order  to  analyze  and 
interpret  the  spectra  a brick  layer  model  (Fig.  la)  is 
usually  applied  assuming  cubic  grains  and  laterally 
homogeneous  as  well  as  identical  grain  boundaries. 
However,  these  brick  layer  assumptions  are  more  or 
less  violated  in  real  ceramics  and  a detailed  in- 
vestigation of  the  impact  of  such  deviations  on  the 
overall  impedance  is  still  missing. 

In  this  contribution  an  overview  of  impedance 
effects  caused  by  some  important  deviations  from  the 
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Fig.  1.  Sketch  of  the  brick  layer  model  (a)  and  important  deviations  in  real  materials  as  considered  in  this  contribution:  (b)  non-cubic  grains 
(here  in  2D);  (c)  laterally  inhomogeneous  grain  boundaries  (here  in  3D);  (d)  grain  boundary  properties  differing  from  grain  to  grain 
according  to  the  probability  distribution  of  the  grain  boundary  conductivity  (here  in  2D).  In  the  brick  layer  approach  (a)  highly  resistive 
grain  boundaries  perpendicular  to  the  electrodes  can  be  neglected  leading  to  the  equivalent  circuit  shown  in  (a). 
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brick  layer  model  is  given,  while  a detailed  analysis 
of  some  special  aspects  is  presented  elsewhere  [8- 
10].  Particularly,  we  consider:  (i)  deviations  from  the 
cubic  grain  shape;  (ii)  the  influence  of  imperfect 
contacts  between  the  grains;  and  (iii)  the  impact  of  a 
distribution  of  different  grain  boundaries  in  the 
crystal  (Fig.  1).  Moreover,  some  conditions  are 
discussed  within  which  the  brick  layer  model  is  a 
useful  tool  to  estimate  electrical  properties  of  grain 
boundaries. 


2.  Model  considerations 

The  potential  distribution  within  a polycrystal  is 
obtained  by  solving  Laplace’s  equation  for  the 
electrical  potential  ( <p ).  The  impedance  in  a sample 
with  bulk  conductivity  <xbulk  and  permittivity  £bulk 
can  be  calculated  by  integrating  the  complex  current 
density,  j = - (a  + ia)e)gmd((p),  along  an  electrode. 
Grain  boundaries  are  assumed  to  exhibit  an  effective 
thickness  wgb  and  an  effective  conductivity  <rgb.  The 
microscopic  reason  for  the  blocking  character  of  the 
grain  boundary  (e.g.  core  effects,  space  charge 
effects  or  a second  phase)  is  not  considered  in  this 
work. 

In  the  case  of  the  two-dimensional  calculations  of 
different  grain  boundary  patterns  (Section  3.1) 
identical  grain  boundary  conductivity  and  thickness 
were  assumed  for  all  grain  boundaries.  In  Section 
3.2,  in  which  the  three-dimensional  simulations  of 
imperfect  grain-to-grain  contacts  are  discussed,  the 
properties  vary  along  a single  grain  boundary:  the 
grain  boundary  conductivity  is  assumed  to  be  zero 
for  an  insulating  grain  boundary  phase  (e.g.  a solid 
phase  or  pores,  Fig.  lc)  and  non-zero  for  the 
contacted  region.  In  Section  3.3  (two-dimensional 
simulations)  a variation  in  the  conductivity  from 
grain  boundary  to  grain  boundary  was  considered, 
using  different  probability  distributions  of  crgb  (cf. 
Fig.  Id). 

The  numerical  solution  was  computed  using  the 
finite  element  software  ‘FLUX-EXPERT’  from 
Simulog  (1  Rue  James  Joule,  78280  Guyancourt 
Cedex,  France).  The  grain  boundaries  may  be  taken 
into  account  via  so-called  ‘interfacial  elements’ 
[9,11].  Further  details  on  such  calculations  are  given 
in  Refs.  [9,12], 


3.  Results  and  discussion 

3.1.  The  influence  of  different  grain  boundary 
patterns 

To  examine  the  influence  of  the  microstructure, 
FE-calculations  in  two  dimensions  were  performed. 
A general  principle  may  improve  the  experimental- 
ists’ prediction  of  the  possible  impact  a given 
microstructure  will  have  on  the  impedance:  if  pos- 
sible, current  lines  make  detours  around  hindrances, 
i.e.  resistive  grain  boundaries.  This  leads  to  an 
interesting  observation.  Detours  are  ‘easy’  for  high 
bulk  conductivity,  and  excess  grain  boundaries  are 
avoided.  However,  detours  are  ‘difficult’  for  low 
bulk  conductivity,  and  additional  grain  boundaries 
must  be  passed,  leading  to  a greater  grain  boundary 
resistance  than  for  the  case  of  high  <xbulk  [9].  Hence, 
the  low-frequency  semicircle  of  the  impedance  spec- 
trum, which  is  usually  exclusively  interpreted  in 
terms  of  grain  boundary  properties,  can  also  be 
influenced  by  properties  of  the  grain  interior. 

Such  detour-effects  can  play  an  important  role  in 
real  ceramics,  particularly  if  broad  grain  size  dis- 
tributions exist  or  if  the  grain  size  distribution  is 
spatially  inhomogeneous.  This  is  illustrated  for  a 
sample  with  inhomogeneous  grain  size  distribution, 
as  shown  in  Fig.  2 (consisting  of  some  7-p.m  large 
grains  and  agglomerations  of  smaller  grains  with 
1-p.m  grain  size).  A simple  calculation  of  an  average 
grain  size  by  counting  the  grains  per  area  leads  to  a 
value  of  1.94  |xm.  The  low-frequency  semicircle  as 
obtained  by  finite  element  calculations  is  considera- 
bly smaller  than  predicted  by  the  brick  layer  model 
for  a grain  size  of  1.94  pm  (see  Fig.  2).  As  a 
consequence,  the  calculation  of  a meaningful  average 
grain  size  requires  a weighting  in  favour  of  large 
grains. 

This  can  easily  be  understood  from  a plot  of  the 
magnitude  of  the  d.c.  current  density  (Fig.  3a):  the 
current  by-passes  the  regions  of  high  grain  boundary 
density  and  thus  the  inhomogeneous  current  dis- 
tribution prefers  large  grains.  However,  this  changes 
for  relatively  low  bulk  conductivity  since  detours 
become  more  difficult  and  a much  more  homoge- 
neous current  distribution  results  as  shown  in  Fig. 
3b.  Hence  the  ‘grain  boundary  resistance’  as  ob- 
tained from  the  low-frequency  semicircle,  depends 
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Fig.  2.  Artificial  microstructure  representing  the  case  of  an 
inhomogeneous  grain  size  distribution  (agglomeration  of  smaller 
grains)  and  the  resulting  impedance  spectrum  as  calculated  by  the 
finite  element  method.  For  comparison  the  spectrum  calculated  by 
assuming  a brick  layer  model  of  1.94  ^m  grain  size  is  plotted  as 
well.  The  grain  boundary  properties  are:  trfh  = 4X  10~“  1/fl,  grain 
boundary  thickness  = 2.8  nm,  permittivity  £,st,  = £,l„1u,  bulk  con- 
ductivity <Tbl]|k  = 10’"4  1/fl.  (Please  note,  that  in  2D  the  unit  of  a 
conductivity  is  1/fl.) 

not  only  on  the  conductivity  of  the  highly  resistive 
grain  boundaries  but  also  on  the  bulk  conductivity 
(Fig.  3c).  As  a consequence,  the  temperature-depen- 
dence of  the  low-frequency  semicircle  ( ‘grain  bound- 
ary semicircle’)  yields  an  apparent  activation  energy 
which  can  differ  slightly  from  the  true  grain  bound- 
ary activation  energy.  We  give  a numerical  example 
for  the  case  considered  in  Figs.  2 and  3:  (i)  if 
K c, bulk  =0-9  eV  and  = 1.8  eV  there  is  a 
‘medium’  temperature  regime  with  an  apparent  grain 
boundary  activation  energy  of  1.7  eV.  (‘Medium’ 
temperature  regime  means  temperatures  for  which 
the  sizes  of  the  two  semicircles  are  of  the  same  order 
of  magnitude.)  The  relative  deviation  from  b can 
be  shown  to  be  approximately  proportional  to  ( 1 — 
^c.buik/^c  gb)  and  thus  can  be  rather  large  for  small 
p 

•^act.gb* 


Moreover,  the  low  frequency  semicircle  is  not 
ideal  but  asymmetrically  distorted.  A fit  using  a 
resistor  in  parallel  with  a constant  phase  element 
with  impedance  Q \ito)  (<2,n  = fit  parameters, 
w = angular  frequency)  yields  a value  of  n =0.916 
for  the  left  half  of  the  low-frequency  semicircle 
while  the  right  half  results  in  n=0.99.  It  has  to  be 
emphasized  that  the  depression  is  not  due  to  a 
distribution  of  relaxation  times  (only  one  grain 
boundary  relaxation  time  exists)  but  is  instead  a 
consequence  of  the  frequency-dependence  of  the 
current  lines:  for  low  frequencies  the  current  detours 
the  excess  grain  boundaries,  while  for  high  fre- 
quencies the  interfaces  are  nearly  dielectrically  short- 
circuited  and  a homogeneous  current  distribution 
results.  This  frequency-dependent  switching  of  the 
current  lines  leads  to  modifications  in  the  spectra 
shape  and  has  no  counterpart  in  quasi-one-dimen- 
sional  systems  as  the  simple  brick  layer  model.  The 
high-frequency  semicircle,  on  the  other  hand,  is  not 
affected  by  the  microstructural  pattern  and  reveals 
the  true  bulk  properties. 

On  the  other  hand,  in  many  cases  the  brick  layer 
model  does  allow  reasonable  estimates  of  the  grain 
boundary  properties.  This  is  discussed  in  more  detail 
in  Ref.  [9].  From  the  calculations  therein  (e.g.  a 
pattern  as  shown  in  Fig.  lb)  and  other  calculations 
for  two-dimensional  grain  boundary  patterns  we 
expect  that  the  brick  layer  model  is  a reasonable 
approximation  for  homogeneous  and  relatively  nar- 
row distributions  of  grain  sizes  and  grain  shapes. 
Serious  problems  are  expected,  for  example,  for  a 
very  broad  grain  size  distribution  with  very  large  and 
very  small  grains  or  a spatially  inhomogeneous  grain 
size  distribution  as  sketched  in  Fig.  3. 

3.2.  The  influence  of  imperfect  contacts 

Grain  boundaries  are  often  laterally  inhomoge- 
neous (i.e.  ‘imperfect’)  in  the  sense  that  ideally 
conducting  as  well  as  totally  insulating  interface 
regions  are  present.  To  give  two  examples:  (i)  solid 
grain  boundary  phases  frequently  only  partially  wet 
the  grains  establishing  a diminished  grain-to-grain 
contact,  and  (ii)  nanopores  along  grain  boundaries 
can  cause  an  imperfect  contact  between  two  grains. 
In  such  cases  an  insulating  layer  partly  separates 
neighboring  grains  (Fig.  lc)  and  the  d.c.  current  is 


Fig.  3.  Top:  magnitude  of  the  d.c.  current  density  within  the  two-dimensional  polycrystal  of  Fig.  2 for  two  different  bulk  conductivities:  (a) 
(7t,u,k  = 2X10-4  1/ft;  (b)  cr  , =1(T3  1/ft.  For  other  parameters  see  Fig.  2.  Dark  indicates  high  current  density  and  light  low  current 


density.  Bottom:  resistance  from  the  low-frequency  semicircle  R2  as 

constricted  close  to  the  conducting  grain-to-grain 
contacts.  Consequently  the  d.c.  bulk  resistance  is 
higher  than  the  ideal  bulk  resistance.  Such  imperfect 
contacts  give  rise  to  a low-frequency  semicircle  in 
the  complex  impedance  plane  (Fig.  4)  which  can  be 


a function  of  bulk  conductivity. 

interpreted  in  terms  of  a transition  of  the  bulk 
resistance  from  the  d.c.  value  to  the  a.c.  (ideal)  value 
due  to  the  dielectric  ‘opening’  of  the  insulating  grain 
boundary  phase  capacitor  [8].  Hence  it  is  not  the 
sequence  of  two  serial  regions  of  different  conduc- 
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Fig.  4.  Impedance  spectra  for  a three-dimensional  polycrystal  (bulk  conductivity  trbulk  =10  1 /fi  cm)  with  an  electrode  area  of  1 cm2  and  a 

thickness  of  5 mm.  All  grain  boundaries  are  identical  but  laterally  inhomogeneous,  consisting  of  an  insulating  phase  (thickness  = 2 nm, 
permittivity  = ebulk)  and  different  amounts  ( n ) of  circular  perfectly  conducting  spots  (cf.  Fig.  lc).  The  insets  sketch  cross-sections  of  the  grain 
boundaries  between  two  grains  (cf.  Fig.  lc)  with  contacted  (white)  and  insulating  (black)  regions.  In  all  three  cases  the  fraction  of  contacted 
area  per  grain  boundary  is  kept  constant  at  12.6%.  The  dc  resistance  (11.4  mft)  is  a pure  bulk  resistance  and  is  larger  than  the  ideal  bulk 
resistance  (5  Mfl)  because  of  the  current  constriction  close  to  the  contacted  spots. 


tivities  but  the  frequency-dependent  shift  of  the 
current  lines  that  leads  to  the  low-frequency  semicir- 
cle. Fig.  4 shows  the  resulting  impedance  spectra  for 
different  grain  boundary  contact  geometries,  all  of 
them  exhibiting  the  same  area  fraction  of  insulating 
grain  boundary  phase  (— acontact)  but  different  num- 
bers of  contact  spots  carrying  the  d.c.  current.  The 
so-called  ‘grain  boundary  resistance’  R2  (diameter  of 
the  low-frequency  arc)  strongly  depends  on  the 
number  of  contacts  while  the  high-frequency 
semicircle  is  not  influenced  by  grain  boundary 
properties  and  still  reveals  the  ‘ideal’  bulk  values  as 
expected  for  a single  crystal.  It  has  to  be  emphasized 
that  the  additional  resistance  is  not  due  to  the  contact 
‘bridge’  which  connects  two  grains  (cf.  Fig.  lc)  but 
is  caused  by  the  current  constriction  in  the  grain  bulk 
close  to  the  conducting  grain-to-grain  contacts. 

A quantitative  estimate  of  the  magnitude  of  the 
‘grain  boundary  resistance’  R2  is  possible  according 
to  Ref.  [8] 

R2~— * (1) 

°bulk  V4acontactn  .sample 

In  Eq.  (1)  n is  the  number  of  grain-to-grain 


contacts  per  grain  boundary  and  acont.lct  the  fraction 
of  established  contact  area  (contacted  area/complete 
area),  rrbulk  denotes  the  bulk  conductivity  and  Ls  lc, 
A.  sample  are  Ae  sample  thickness  and  area,  respective- 
ly. Thus  R2  is  determined  by  two  independent 
parameters  (acontact,«),  and  an  evaluation  of  the 
contacted  area  (or  of  the  fraction  of  ‘blocked’ 
current)  solely  from  the  impedance  spectrum  is  not 
possible.  The  ‘grain  boundary  capacitance’,  on  the 
other  hand,  allows  an  estimate  of  the  thickness  of  the 
insulating  grain  boundary  phase  [8], 

To  help  estimate  to  what  extent  current  constric- 
tion might  cause  grain  boundary  impedances  in  real 
electroceramic  materials,  a plot  is  given  (Fig.  5) 
which  relates  i?2  //? , (i.e.  the  ratio  of  the  semicircle 
diameters)  to  the  fraction  of  contacted  area  (acon,.lct) 
for  varying  contact  numbers  per  grain  boundary. 
From  Fig.  5 and  also  from  the  estimate  according  to 
Eq.  (1),  it  can  be  seen  that  pronounced  low-fre- 
quency semicircles  require  rather  small  contact  areas, 
particularly  if  many  contact  spots  are  involved.  In 
many  cases  we  expect  the  grain  boundary  resistance 
due  to  current  constriction  to  be  of  the  order  of 
magnitude  of  the  ideal  bulk  resistance  or  even 
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Fig.  5.  Ratio  of  R2IR]  (‘grain  boundary  resistance’ /ideal  bulk 
resistance)  depending  on  the  fraction  of  contacted  area  aconlact,  for 
different  numbers  of  contacts  per  grain  boundary  n.  The  plots  are 
calculated  using  a more  generally  valid  estimate  than  Eq.  (1)  as 
given  in  Ref.  [8]. 


smaller.  However,  in  porous  material  (e.g.  freshly 
pressed  samples)  or  if  there  are  only  very  few  gaps 
(pin-holes)  in  a highly  resistive  grain  boundary  phase 
current  constriction  effects  can  be  considerably 
larger. 

Indications  that  current  constriction  effects  are 
relevant  can  be  obtained  from  the  temperature, 
partial  pressure,  bias  and  grain  size  dependence  of 
R2.  If  current  constriction  plays  a role,  the  tempera- 
ture and  partial  pressure  dependence  of  R2  should  be 
close  to  that  of  the  bulk  as  already  predicted  by 
intuitive  models  [3,5-7,13],  A bias-dependent  ‘grain 
boundary  resistance’,  on  the  other  hand,  suggests 
other  mechanisms  than  current  constriction.  The 
grain  size  dependence  of  the  grain  boundary  resist- 
ance, finally,  is  expected  to  deviate  from  the  inverse 
linear  relationship  suggested  by  the  brick  layer 
model.  Further  information  is  given  in  Refs.  [8,12], 

However,  the  situation  becomes  more  complicated 
if  partially  blocking  and  totally  blocking  interface 
regions  exist  simultaneously.  For  example,  an  in- 
sulating grain  boundary  phase  partially  wets  the 
grains  while  at  the  established  grain-to-grain  contacts 
an  additional  transfer  resistivity  exists.  Current  con- 
striction effects  and  ‘true’  grain  boundary  contribu- 
tions due  to  the  transfer  process  superimpose  and 
lead  to  a rather  complex  situation. 


To  get  a first  insight  into  this  topic  we  consider  the 
following  grain  boundary  case.  A solid  grain  bound- 
ary phase  (thickness  2 nm)  partially  wets  the  grains, 
leaving  a circular  grain-to-grain  contact  with  a 
diameter  of  54%  of  the  grain  boundary  length  (2 
p,m)  (as  sketched  in  Fig.  lc).  At  the  established 
contact  we  assume  a transfer  hindrance  which  can  be 
described  by  a 2-nm  thick  ‘layer’  of  conductivity  <rgb 
and  permittivity  £bulk.  As  shown  in  Fig.  6,  again  one 
low-frequency  semicircle  arises  in  the  complex 
impedance  plane.  However,  this  low-frequency 
semicircle  includes  both  bulk  contributions  (due  to 
current  constriction)  and  ‘true’  grain  boundary  con- 
tributions owing  to  the  partial  blocking  at  the 
established  contact.  Hence  a ‘mixed’  activation 
energy  results  and  even  a change  of  the  activation 
energy  with  temperature  from  the  bulk  value  to  a 
higher  or  lower  value  (non-Arrhenius-behaviour)  is 
possible.  Nevertheless,  the  contribution  of  current 
constriction  within  the  grains  can  be  estimated  by 
Eq.  (1). 

Such  ambiguous  situations  arise  if  the  capacitance 
of  the  insulating  grain  boundary  phase  is  larger  than 
the  capacitance  related  to  the  transfer  process  at  the 
established  contact.  Since  relevant  current  constric- 
tion effects  require  relatively  small  contacted  areas 
(see  Fig.  5)  we  expect  that  this  condition  for  the  two 


♦ transfer  resistance  at  contact 
o ideally  conducting  contact 
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Fig.  6.  Impedance  spectra  for  a three-dimensional  crystal  (bulk 
conductivity  <zbu|k  = 10  7 1/ft  cm)  with  and  without  an  additional 
transfer  resistance  at  the  established  grain-to-grain  contact.  All 
grain  boundaries  are  identical  and  exhibit  one  contact  spot  with 
contact  area  aconucl  = 0.23.  Parameters:  electrode  area,  1 cm2; 
sample  thickness,  5 mm;  thickness  of  the  insulating  phase,  2 nm; 
permittivity  of  the  insulating  phase,  ebulk;  partial  blocking  at  the 
established  contact  according  to  a 2-nm  thick  ‘layer’  of  con- 
ductivity, £rgb;  and  permittivity,  ebu]k. 
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capacitances  is  often  fulfilled  at  real  grain 
boundaries.  However,  if  the  opposite  situation  is 
valid  (capacitance  of  the  insulating  layer  smaller  than 
the  transfer  capacitance)  a third  semicircle  appears, 
and  an  approximate  separation  of  current  constriction 
and  true  grain  boundary  resistance  is  possible,  as 
shown  for  a similar  problem  (fuel  cell  cathode  with 
the  current  being  constricted  close  to  the  three  phase 
boundary  region)  in  Ref.  [14], 

3.3.  The  influence  of  the  distribution  of  different 
grain  boundaries  in  space 

Even  if  each  individual  grain  boundary  is  laterally 
homogeneous,  the  structural  and  chemical  qualities, 
and  thus  the  electrical  properties,  may  vary  from 
boundary  to  boundary  leading  to  a certain  distribu- 
tion of  grain  boundary  conductivities  as  sketched  in 
Fig.  Id.  One  example  where  such  a distribution  is 
expected  is  for  the  case  when  the  grain  boundary 
resistance  is  caused  by  a space  charge  depletion 
layer.  Due  to  an  exponential  relation,  relatively  small 
variations  of  the  space  charge  potential  result  in 
strong  variations  of  the  effective  grain  boundary 
resistivity. 

We  probed  the  influence  of  distributions  of  differ- 
ent grain  boundaries  in  two  dimensions  for  some 
probability  distributions  of  the  grain  boundary  con- 
ductivity crgb  (Figs.  7-9).  Following  the  given 
distribution  the  grain  boundary  conductivities  were 
randomly  distributed  in  the  square  microstructural 
pattern  (cf.  Fig.  Id,  the  calculations  are  performed 
for  a sample  consisting  of  10  X 10  grains). 

Fig.  7 shows  the  spectrum  for  a discrete  Gaussian 
distribution  of  log(cr,b)  calculated  by  the  finite 
element  method  and  a spectrum  according  to  the 
brick  layer  model  using  the  ‘mean’  grain  boundary 
conductivity  (1(T7  1/fl).  Such  a distribution  of 
log(crgb)  corresponds  to  a Gaussian  distribution  of 
the  space  charge  potential.  Even  though  the  extreme 
conductivity  values  of  the  distribution  differ  by  a 
factor  of  100  the  brick  layer  model  allows  a rather 
satisfying  estimate  of  the  ‘mean’  grain  boundary 
properties.  However,  again  an  inhomogeneous  cur- 
rent distribution  results  which  depends  weakly  on  the 
bulk  conductivity.  Therefore,  bulk  properties  slightly 
influence  the  ‘grain  boundary  semicircle’  (cf.  Fig.  3), 
and  a minor  difference  between  the  activation  energy 


Fig.  7.  Calculated  impedance  spectrum  of  a polycrystal  (10X10 
grains)  for  a discrete  Gaussian  distribution  of  the  grain  boundary 
conductivities  (cf.  Fig.  Id).  The  Gaussian  distribution  is  on  a 
logarithmic  scale.  The  corresponding  brick  layer  model  assumes  a 
conductivity  according  to  the  ‘mean’  value  of  the  distribution 
(10  7 1/fl).  Parameters:  two-dimensional  bulk  conductivity,  10  4 
I/O;  grain  size,  2 p.m;  grain  boundary  thickness,  2 nm. 
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Fig.  8.  Calculated  impedance  spectrum  of  a polycrystal  (10X10 
grains)  for  a very  broad  distribution  of  the  grain  boundary 
conductivity  (cf.  Fig.  Id).  The  corresponding  brick  layer  model 
assumes  a conductivity  according  to  the  center  of  the  distribution 
(10  1/fl).  Other  parameters  as  in  Fig.  7. 
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Fig.  9.  Calculated  impedance  spectra  of  a polycrystal  (10X10 
grains)  for  different  bimodal  distributions  of  the  grain  boundary 
conductivity  (cf.  Fig.  Id).  While  70%  of  the  grain  boundaries 
exhibit  bulk  conductivity  trbulk  = crgb  , = 10-4  1/ft,  varies 

between  10~7  1/ft  and  zero.  Other  parameters  as  in  Fig.  7. 


of  the  low-frequency  semicircle  and  the  true  activa- 
tion energy  of  the  grain  boundary  resistivity  can 
occur.  Moreover,  the  low-frequency  semicircle  ex- 
hibits a weak  asymmetric  depression  (Fig.  7). 

The  situation  becomes  more  difficult  for  a very 
broad  conductivity  distribution  of  constant  probabili- 
ty (Fig.  8).  A strongly  deformed  ‘grain  boundary 
semicircle’  results  and  a brick  layer  model  using  the 
conductivity  of  the  center  of  the  distribution  (10~7 
1/fi)  underestimates  the  ‘grain  boundary  resistance’ 
considerably.  On  the  other  hand,  a serial  connection 
of  grain  boundaries  according  to  the  given  probabili- 
ty distribution  yields  a much  larger  value  of  R2  than 
the  two-dimensional  distribution  (222  kfl).  This 
again  demonstrates  the  influence  of  detours  around 
high-resistive  grain  boundaries  to  lower  the  resist- 
ance but  also  the  problems  to  deduce  a meaningful 
parameter  from  the  ‘grain  boundary  resistance’. 
Furthermore,  the  ‘grain  boundary  resistance’  is  dis- 
tinctly influenced  by  the  bulk  conductivity.  As  in  the 
case  of  different  grain  boundary  patterns  (see  Section 
3.1),  detour  effects  around  very  blocking  grain 
boundaries  contribute  to  the  low-frequency  ‘semicir- 
cle’. Consequently  the  apparent  activation  energy  no 


longer  reflects  the  activation  energy  of  the  grain 
boundary  conductivity  even  if  all  crgb  of  the  dis- 
tribution have  identical  activation  energies. 

Let  us  finally  discuss  the  extreme  case  in  which 
70%  of  the  grain  boundaries  are  ideal  and  exhibit 
bulk  conductivity  while  30%  are  highly  resistive 
(bimodal  distribution).  Fig.  9 demonstrates,  how  the 
‘grain  boundary  arc’  depends  on  the  conductivity  of 
the  highly  resistive  interfaces  (crgb  2).  For  moderate 
grain  boundary  conductivity  the  arc  is  semicircle-like 
and  depends  mainly  on  a b 2 but  only  slightly  on  the 
bulk  conductivity.  For  highly  blocking  grain 
boundaries,  on  the  other  hand,  the  low-frequency 
‘semicircle’  becomes  distinctly  distorted  and  its  size 
depends  only  weakly  on  crgb  2 but  mainly  on  the  bulk 
conductivity.  In  the  extreme  case  of  30%  totally 
insulating  grain  boundaries  the  low-frequency  arc 
reaches  a maximum  and  the  so-called  ‘grain  bound- 
ary resistance’  turns  out  to  be  a part  of  the  d.c.  bulk 
resistance  caused  by  detours  around  the  very  block- 
ing grain  boundaries.  Hence  the  mechanism  deter- 
mining the  resistance  of  the  low-frequency  arc 
changes  with  decreasing  grain  boundary  conductivity 
and  the  apparent  activation  energy  varies  from  the 
activation  energy  of  crgb  2 to  the  bulk  value  for 
decreasing  cr£b  2/crbulk  ratio.  This  is  demonstrated  in 
Fig.  10  for  bulk  and  grain  boundary  activation 
energies  of  0.6  and  1.2  eV,  respectively. 

These  examples  demonstrate  that,  on  one  hand,  a 
brick  layer  analysis  can  yield  meaningful  values  as 
long  as  relatively  narrow  Gauss  distributions  of  grain 
boundary  properties  are  considered,  while  it  can 
completely  fail  in  other  cases.  A more  detailed 
investigation  in  three  dimensions  is  in  progress  [10] 
and  will  reveal  more  quantitatively  the  conditions 
under  which  a brick  layer  analysis  is  useful. 

Finally  we  can  conclude  from  Sections  3.1  to  3.3 
that  two-  or  three-dimensional  current  lines  can 
depend  on  bulk  conductivity  and  a.c.  frequency.  This 
leads  to  new  features  which  do  not  exist  in  related 
one-dimensional  problems  as  for  example  to  bulk- 
dependent  ‘grain  boundary  semicircles’  and  to  dis- 
torted or  additional  arcs  only  due  to  a change  of  the 
current  lines  with  increasing  frequency.  As  a conse- 
quence, the  resistance  of  the  low-frequency  semicir- 
cle ( R2 ) does  frequently  not  represent  the  exact 
resistance  of  the  grain  boundaries  themselves.  On  the 
other  hand,  the  term  ‘grain  boundary  resistance’  is 
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Fig.  10.  Bulk  (/?,)  and  ‘grain  boundary  resistance’  ( R: ) for  the 
bimodal  distribution  of  grain  boundary  properties  as  shown  in  Fig. 
9.  The  bulk  conductivity  tr1>ul).  has  an  activation  energy  of  0.6  eV 
while  cr,h  , is  activated  with  1.2  eV.  Notice  that  R)  '(T)  reflects  the 
activation  energy  of  trfK2  at  high  temperatures,  but  R,'(T) 
approaches  bulk  activation  at  low  temperatures. 


still  meaningful  in  the  sense  that  R2  is  caused  by  the 
grain  boundaries  and  does  not  exist  in  a corre- 
sponding single  crystal. 

However,  in  all  cases  the  high-frequency  semicir- 
cle reveals  the  ideal  bulk  values  as  expected  for  a 
single  crystal,  since  for  the  entire  high-frequency  arc 
the  grain  boundaries  are  dielectrically  conducting 
and  hence  negligible.  This  is  even  true  if  the  d.c. 
resistance  is  a pure  bulk  resistance,  i.e.  if  the  low- 
frequency  semicircle  is  caused  by  detours  around 
blocking  grain  boundaries. 

4.  Conclusions 

• Deviations  from  the  cubic  grain  size  are  accept- 
able for  a brick  layer  analysis  as  long  as  the  grain 
size  distribution  is  relatively  narrow  and  spatially 
homogeneously,  and  as  long  as  no  pronounced 
anisotropy  in  the  grain  shape  distribution  exists. 
However,  agglomerates  of  small  grains  or  very 


broad  distributions  can  lead  to  considerable  er- 
rors. 

• If  imperfect  contacts  (lateral  inhomogeneities) 
between  grains  occur,  current  constriction  in  the 
grain  becomes  important.  This  leads  to  an  appar- 
ent grain  boundary  semicircle,  the  resistance  of 
which  is  purely  bulk-dependent.  An  estimate  of 
the  additional  resistance  can  be  given  according 
to  Eq.  (1),  however,  a quantitative  determination 
of  the  fraction  of  blocked  grain  boundary  area 
solely  from  the  impedance  spectrum  is  not  pos- 
sible. A simple  brick  layer  analysis  fails  in  such 
cases. 

• The  brick  layer  model  can  be  used  to  estimate 
mean  grain  boundary  properties  if  a relatively 
narrow  (Gaussian)  distribution  of  grain  boundary 
conductivities  is  given.  However,  very  broad  or 
bimodal  distributions  lead  to  large  deviations  and 
a meaningful  analysis  using  the  brick  layer  model 
fails.  In  such  cases  the  low-frequency  arc  can  be 
distinctly  influenced  by  bulk  properties  and  con- 
siderably distorted.  This  can  lead  to  apparent 
activation  energies  which  differ  from  the  true 
grain  boundary  activation  energy. 

• Grain  boundary  arcs  deviating  from  ideal  semicir- 
cles do  not  only  occur  for  distributions  of  relaxa- 
tion times,  but  also  for  microstructures  which 
differ  from  the  brick  layer  model  even  if  only  a 
single  grain  boundary  relaxation  time  is  involved. 

• In  all  cases  considered  the  high-frequency 
semicircle  still  reveals  the  ideal  bulk  resistance 
and  capacitance  as  expected  for  a single  crystal. 
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